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ON THE CIRCULATION MANIFOLD FOR TWO ADJACENT LIFTING SECTIONS 

LUCA ZANNETTU AND ANGELO IOLLO* 


Abstract. The circulation functional relative to two adjacent lifting sections is studied for two cases. In 
the first case we consider two adjacent circles. The circulation is computed as a function of the displacement 
of the secondary circle along the axis joining the two centers and of the angle of attack of the secondary 
circle. The gradient of such functional is computed by deriving a set of elliptic functions with respect both 
to their argument and to their period. In the second case studied, we considered a wing-flap configuration. 
The circulation is computed by some implicit mappings, whose differentials with respect to the variation of 
the geometrical configuration in the physical space are found by divided differences. Configurations giving 
rise to local maxima and minima in the circulation manifold arc presented. 

Key words, potential flow, elliptic functions, optimal shape design. 

Subject classification. Fluid Mechanics, Applied and Numerical Mathematics 

1. Introduction. Considering the growth of research and industrial interest in aerodynamic shape 
optimization, we felt that a reference test case amenable to analytic solution and yet complex enough to be 
representative of problems arising from real world applications, was lacking. For this reason we propose in 
this paper two optimization test cases based on the potential flow solution. The first is represented by the 
problem of finding an extremum of the circulation manifold relative to the flow about two circles, when the 
secondary circle and its trailing edge are displaced. The global circulation and its gradient arc computed 
by analytical tools. In the second test case, the problem of computing the circulation and its gradient for 
the flow about two airfoils, when the second airfoil is translated and rotated with respect to the first one, is 
solved. 

The solution is determined by a classical analytical technique that requires two steps. The first step, 
which is common to both test cases, is the determination of the flow solution in the plane of the two circles. 
In the plane of the airfoils only 3 parameters arc varied: the position of the leading edge of the flap (X S ,Y S ) 
and its deflection 0. On the plane of the circles it is shown that there are 6 parameters changing as the flap 
is moving in the physical space. The six parameters are the modulus and direction of the speed at infinity, 
the radii of the the two circles and the angles defining the position of the two trailing edges. Using the 
Lagally [1] solution for the flow about two circles, it was possible to compute the analytical differential of the 
total circulation about the two circles. The second solution step is the transformation of the two airfoils into 
two circles using a Thcodorsen transform for the primary airfoil and a Garrick transform for the other. Any 
two airfoils can be transformed by this technique. Yet, the mappings used are known implicitly by means of 
a FFT. Therefore, for simplicity, the Jacobian of these transformations with respect to the translation and 
rotation of the flap is computed by divided differences. Such differentiation can be conducted to any order 
of accuracy at a very limited computational cost. 
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Projet SINUS, Sophia Antipolis Cedex, France. 
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Optimization consists in defining a functional which represent the cost (or merit) of a solution, as 
a function of the geometrical configuration. The optimal solution is the minimum (or the maximum) of 
such a functional. An optimization algorithm is based on three ingredients: the flow solver, a routine to 
compute the gradient, and an optimizer, i.e., a strategy to march toward the extremum of the functional, 
sec for example [2] and references therein. Using the exact computation of the gradient that we propose 
it is possible to validate the accuracy of numerical gradients obtained by adjoint or sensitivity analysis 
formulations. Furthermore, general type descent algorithms can be tested using the exact gradient given in 
this paper, in order to evaluate performance in the context of aerodynamic optimization. For example, the 
effects of local minima on the performance of various optimizers can be investigated. There was another 
attempt at computing an exact test case for the circulation about wing-flap section configurations [3]. It 
was found that the derivation presented therein is based on some erroneous assumptions. In fact, because 
of the motion of the flap, the primary airfoil shape and the flap shape are deformed, making that test case 
of no practical relevance. 

Section 2 is devoted to the study of the analytical gradient of the circulation about two adjacent circles 
with given trailing edge. A first exact test case is presented. In Section 3 the mapping from the two airfoils 
plane to the two circles plane is sketched. Finally, he results relative to the second test case proposed, 
obtained by divided differences, arc presented. 

2. The circulation manifold of two adjacent circles. In this section we will assume that the 
two profiles arc indeed two circles for the reason that in this case it is possible to determine by analytical 
differentiation the dependence of the circulation and hence of the lift, as a function of some geometrical 
parameters. 

Consider two circles (Fig. 2. 1-2.3), the first is the unit circle centerd at the origin, the second has radius 
r s and center at x = x s on the real axis of the plane z — x + iy. The two circles arc immersed in a potential 
incompressible flow. The speed at infinity of such flow is 9oo = |^|e ia ° and the Kutta condition is imposed 
at the points z v and which belong to the primary and secondary circle respectively. They arc determined 
by ip p — arg( 2 p ) and <p s = arg(z 5 ). The solution to this flow field is due to Lagally [1] and Ferrari [4]. 



Let 

( 2 . 1 ) 


Fig. 2.1. Plane of the two circles. 


b = 


r ; + 1 


2x s 


c = \/ (x s - b) 2 - rl 
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Fig. 2.3. Plane of the rectangle B. 

The two circles considered belong to a family of Apollonius circles with focal points at c + b and c — b 
on the real axis. The mapping 


transforms this Apollonius manifold to another, whose focal points arc at v = — 1 and v = 1. 

Let us exploit the invariance of the complex potential w under conformal mappings, and let us consider 
the flow on the transformed plane v. The complex velocity at infinity is 


Qlc = \Qc 


dw dw dz , , •_ 

e 00 = lim — = lim — — = ck/^e 0 
v-+oc dv z —* oo dz dv 


The bilinear transform 


maps the two circles in the plane 2 onto two concentric circles with center in the origin. The internal circle, 
whose radius is r* < 1, corresponds the primary circle, the external, whose radius is r e > 1, corresponds to 
the secondary circle. The region bounded by the two circles corresponds to the region where the flow field 
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takes place on 2 and the point ji = 1 corresponds to infinity on the plane 2 as well as on v. The trailing 
edges Zp and z s are mapped to // p = r* cxp(h? p ) and fi s = r e cxp(ii3 s ) respectively. 

Now, the mapping 

(2.5) A = log !i 


transforms the region bounded by the concentric circles into a periodic row of rectangles in the plane 
A = £ + ir]. We denote by D the rectangle —(3 < £ < a and — 7T < 77 < 7 r, with (3 ~ — logr^ and a = logr e . 
The period of this row of rectangles is 2 z 7T. The primary circle goes onto the segment f = — /?, — n < 77 < rr, 
while the secondary circle goes onto the segment £ = a, — 7 r < r] < 7r. The trailing edges are mapped to 
A p = —j3 -f ifip , and A s = a -f 

The complex speed inside B is singular only at the origin A = - 0, which corresponds to infinity in v and 
2 . The possible singularities can be a pole of order 2, responsible for the far field speed on the physical 
plane, and a pole of order 1, relative to a possible vortex at infinity whose intensity is equal to the opposite 
of the circulation about the two circles. In addition, the complex speed, which is periodic in the direction of 
the imaginary axis with period 2z7r, must be periodic in the direction of the real axis with period 2 (a + /3) 
in order to satisfy the impermeability condition at the walls of the two circles. As a matter of fact, the 
impermeability condition is obtained by successive reflections with respect to the two sides of the rectangle, 
so that it is periodic with period 2 (a + (3) in that direction. Therefore, the period is formed by B and the 
mirror image of B. We conclude that the complex speed is expressed by an elliptic function with semi-periods 
lj — a + 3 and u> = Z7r, and that it has poles of order 1 and 2 located in A = 0 and A = —2(3. Hence, the 
complex speed is 

d\n iV 

(2-6) — = - — I [C(A) - C(A + 2/?)] - 2[Q>(A) - Q^A + 2/3)] + in 

and, therefore, the complex potential is 

(2-7) w = ~ log + 2[Q^C(A) - Q, C(A + 2/3)] + in A 

where £, p, <r , arc Weicrstrass functions with semi-periods uj = a + (3, Q = in. 

The constants T and k are respectively the global circulation and the circulation about the secondary 
circle. They are determined by imposing the Kutta conditions 


( 2 . 8 ) 

We have 

(2.9) 

where 



= 0 


— liri(C s ~ C p ) 
Ap A s 


K — 



- 2 iC p 


C s =Q^p(A s )-Q ooP (A s - 2/3) 
C p — QocP(^p) — < ?c»p(Ap + 2/3) 
A s = C(A S ) - C(A S + 2/3) 

(2.10) A p = C(A P ) — C(A P + 2/3) 
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2.1. Differentiation. Two adjacent airfoils can always be mapped using the Theodorsen- Garrick trans- 
form onto two circles belonging to an Apollonius family of circles on the complex plane i/, with focal points 
at v = —1 and v = 1. When either the geometric or the flow configuration of the two airfoils changes, 6 
parameters may change on the plane of the two circles 


( 2 . 11 ) 


\Qoo\, A), a, /3, tip 


These parameters, in turn, determine the circulation T (see Eq. (2.9)) and hence the lift. Consider the 
problem of two adjacent rigid airfoils, the lift is a function of 4 parameters: the angle of attack and the three 
degrees of freedom of the second airfoil (every geometric or flow configuration can be reduced to the case of 
fixed primary profile and = 1). Therefore keeping constant the angle of attack at oo, in the easiest case 
of two circular airfoils the degrees of freedom to be considered are 3: the distance between the centers x s , 
and the position of the two trailing edges <p p . In such simple case the circulation (and the lift) is a 3 
parameter manifold in a 6-dimensional space represented by the parameters (2.11). Therefore, in general we 
can express all derivatives of the circulation as a combination of the derivatives of Eqs. (2.10) with respect 
to the parameters (2.11). 

Since u> = a + /3, to differentiate with respect to a or /3 implies differentiation with respect to the 
semiperiod of the elliptic functions oo. These derivatives are given in Section 2.2. 


2.1.1. Derivative d Q 

dC , 


da 


= q: 


P^As) + 


dp{ K) 


-Q c 


dp(X P ) 

da Wo ° dco 


duo 

dp(X p + 2/3) 


n\ - n a\ . &P(^s + 2/3) 
P (X s + 2/3) -I — 


duo 


dA s 

da 


‘p(A$) "f 


ogxs) 

duo 


— p(A 5 + 2/3) + 


dC(A$ + 2/3) 
duo 


( 2 . 12 ) 

2 . 1 . 2 . 


dAp = dgx p ) _ dgx p + 2 / 3 ) 

da duo duo 


Derivative dp. 


9£ 1 = n . n 

d(3 du V ° 


2p'(\ s + 2/3) + W »+ 2/?) 


(2.13) 

2.1.3. 


L.'a 

QQ ~ Qoo P (”p) "3" 

dA s = ac(A s ) 
dp dui 

dAp = p(A p ) + 9C(Ap) 


dp(A P ) 

duo 

2p(A s + 2/3) + 


p '( Ap + 2/ 3) + M^±M' 


Q 

d((\ s +2p) 


duo 


dp 

Derivative d # 


duo 


— p(A p + 2/3) -f 


^C(A P + 2/3) 
duo 


dC 

of = i\Q*ooP'( a.) - QooP’iK + 2/3)1 

dd 


dd 

dA 


= 0 


(2.14) 


~ *[ — p(A s ) + p(A s + 2/3)] 

dA v 


di) 


^ = 0 
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2.1.4. Derivative d$ p . 


2.1.5. Derivative du 


(2.16) 


2.1.6. Derivative dn 


i{Q' <x p , (\p)-Q ooP '{\p + 2f3)} 


— i[— p(Ap) + p(Ap -f 2 ft)] 


~ _ — ^[Q^op(As) + Qoop{K + 2/3)] 
o~ = — ^[QooP(^p) d- Qoop(A p + 2/3)] 


(2 - 17) ^ = ° 

2.2. Derivative of the Weierstrass functions with respect to the period. When the geometric 
configuration varies on the physical plane, the radii r,, r e change as well. Therefore on the plane A the semi- 
period u) = a 4* /3 of the elliptic functions p(u; u>,u>) and ((u;u;,d’) of Weierstrass varies when the geometric 
configuration is changing. It follows that in the derivatives of Eqs. (2.12-2.13) with respect to a and ft, 
besides the partial derivatives of p(u;u, u>) and £(u;u;,u)) with respect to u , there appear the derivatives 
with respect to lj, For example, consider the derivative with respect to a of p(A s ), appearing in d a C s . Since 
X s = a + we have 

dp(u = a + ii9 s ; cu — a 4- /3, u) = in) dp(u;cv,u>) du dp(u\u,cj) du) 

da du da dw da 

We mention in passing that 


(2.18) 


a C(u) , ^ dp(u) „ N 

— — = -p(u), -x— =PW 


du du 

On the contrary, the derivation of these functions with respect to the semi-period is not as easy. 

It is convenient both for the computation and for the derivation to express £(u) and p(u) in terms of 
Jacobi ?3 a (u) functions [5]: 


(2.19) 




p(u) = ei + 


1 K(0) ih(v) 

4 u) 2 ^2(0)i?i'v) 
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with 


2c o 

1 0J"(O) 
V ~ I2w&' 1 ( 0 ) 

e.4n 2 W 


(2-23) = [^(0)+^(0)] 

where denotes derivation with respect to v. Jacobi (t>) functions can be expressed as 

oo 

Mv) = 2g< ^(-l)"g"( n+1 ) sin[(2 n + 1)™] 


tf 2 (iO = 2g* < ?" ( " +1) cos[(2n + 1)™ 

71=0 

OO 

^3(v) = 1 + 2 ^ <7" cos(2n7rt?) 

n=l 

oo 

^ 4 (t;) = 1 + 2 ^^(— l) n q n cos(2n7n?) 


id 

r = — 

Cl? 


The convergence rate of the series in Eqs. (2.24) is high when |cj| < |cl?| and ) >0. It is always possible 
to satisfy such conditions with a proper choice of the rectangle periods. Up to now we always assumed 
u — a + (3, u — in, and hence the first condition is not always met. If it turns out that |ct>| > 1^1, it is 
possible to use the allowed transformations for the periods of elliptic functions (sec [5]), as the periods are 
not uniquely defined. In our case we can define new semi- periods id and u? which are obtained from the 
previous by means of a transformation of first order 


so that the conditions mentioned are both satisfied. 
Finally, we obtain 

d((u) = 1 

did 2Cl^i(v) 


L j m , ggiM 

\{v)[ dw 


mv ) 

did 


1 ^i(u) \ u dr] 

2C1 2 $i(n) + cl? 2 / ^ duj 


dp(u) de 


= ~ + 2(p{u) - ei) -- + 


did did 


+2 {p(u) - ei) 


J dd\{ 0) 1 <9tf 2 (0) 1 


a du $[(0) did tf 2 (o) 






where Cl — id — a ~ (3 and J = 1 if |w| < |u;|, whereas Cl = id = in and J = 0 if |u;| > \<d\. 
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2.3. Differentiation of //(z). Let ji{z) be the chain of trails fori nations which leads from the plane of 
the two adjacent circles z to )i. the plane of the concentric circles. Keeping a 0 and tp p fixed, the mapping 
is uniquely defined as a function of x s and ip s . We derive the parameters (2.11), keeping a fixed velocity at 
infinity, with respect to x s and tp,. These derivatives involve a, 6, r, and r e . It is r. = p(- 1), r e = fi(x 3 +r s ). 


(2.29) 


db = / 2 x, _ J_\ 

dx s Vx^-r^ + l x s ) 


(2.30) 


(2.31) 


(2.32) 


dc = W 1 fa) 

dx s c 

db dc db dc 

da __ dx 3 dx s dx s dx s 

dx s x a + r s - b + c x s + r 8 - b — c 

db dc db dc 

d{3 _ dx s dx s dx s dx s 

dx s c — b — 1 c + 6 + 1 


(2.33) 


(2.34) 


(2.35) 


whereas 

(2.36) 


d0 s _ /3A. 

aZ = ,ra (S 


;) 


7m 


/ db dc 
dx s dx 9 
z s — b + c 


1 - 


06 0c \ 

dx s dX a I 


b — c 


dx. 


0x 5 / 


( db dc db dc \ 

dx s dx, dx s + dx, 


c 6 + z v 


V 


-c — b + z p 


) 


do. 

dipt 

Im 


■ ,m (fe) 


r 3 (— sin<p 3 + i cos <p s ) r 8 (- suup s + 


z s — b + c 


z« — b 


-i cos 


0a _ 0/3 _ d$p _ 

df s ~ dO s ~ d0 s - 


2.4. Circulation manifold. As we change x s and 6 S there are four remarkable situations. The lift may 
asymptoticaly increase or decrease with increasing distance between the two circles. In addition there may 
be a local maximum, or a local minimum. The variation of 6 S dees not give rise to any interesting change in 
the shape of the manifold, except in Fig. 2.10 where a local maximum of the circulation with respect to 6 S is 
shown. The different behaviors arc due to the interplay between the radius of the secondary circle and 0 P , as 
it is shown in the following figures. This is basically the reason for which this test case may be interesting: 
gradient based method would be not able to exit a local minima, whereas genetic algorithms can, see for 
example [7] . The case shown in Figs. 2. 4-2. 5 is for a secondary circle of radius r s = 0.2, 9 P = -10°, ao = 10°, 
1.3 < x s < 6.3, -10° < 0 S < 10°. Each of these intervals is divided in 20 segments and the circulation and 
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its derivatives are evaluated at every resulting grid point. For this configuration lift is increasing with the 
distance, whereas in the case of Figs. 2.6-2. 7 lift is decreasing with distance (r s = 0.6, t 9 P = 30°, a 0 = 10°, 
1-0 < x s < 6.6, —30° < 0 S < 30°). In these figures and in the following the results are presented with respect 
to the grid points rather than explicitly with respect to the variable relative to the axis. 

In Figs. 2. 8-2. 9 it is shown how with a proper choice of the parameters it is possible to obtain a local 
minimum (r s = 0.05, 0° p = -30, ag = 10, 1.15 < < 6.15, -30° < 0 a < 30°). In contrast, with r s = 0.05, 

Op ~ 30, ag = 10, 1.15 < x s < 6.15, —30° < 0 S < 30°; we have a local maximum, with respect to both x s 
and 0 S (see Figs. 2.10- 2.12). 

We limited ourself to the explicit computation of the derivatives with respect to x s and 6 S in order to 
present conveniently the results and to keep a link with real world applications where these are basically 
the parameters to be varied. It is straight forward, using the derivatives of the elliptic functions, that we 
evaluate in the sections before, to determine the differential with respect to any other geometrical quantity 
on the plane of the adjacent circles. 



Fig. 2.4. Circulation as a function of x s (xs) and 0 S (ths): the circulation increases unth x 3 . 


3. The Theodorsen-Garrik transform revisited. Here we adapt the Theodorsen-Garrik transform 
as proposed by Ives [6], in order to differentiate it by divided differences. 

Our scope is to transform the domain exterior to two adjacent airfoils (sec Fig. 3.1) on the plane Zf , 
into the region bounded by two concentric circles on the plane z c . Let us first transform the primary profile 
into a nearly circular region by the Karman-Trefftz transform 



with r = 7r/(27r - <S P ), S p being the trailing edge angle of the primary profile, z Tp is its location, and zn v is 
a point in the vicinity of the leading edge, inside the airfoil. By the Theodorsen transform the quasi circle 
on C is mapped onto a unit circle on z 


(3.2) 


C z ex p{ ^~^ [( a j "i - ibj) z ^1} 

3=0 
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Fig. 2.5. Derivative of the circulation with respect to x 8 . 



Fig. 2.6. Circulation as a function of x 3 and 9 S : the circulation decreases with x 8 - 

The coefficients dj and bj of the suitably truncated series in Eq. (3.2) are iteratively found by an FFT, 
imposing that the points of the unit circle on z are mapped on tbe border of the nearly circular domain on 
For the details see [6]. 

The sequence of mappings that transformed the primary prohle into a unit circle acted on the secondary 
profile as well. Yet, the tangent discontinuity at the trailing edge of the secondary profile is still there, 
therefore it is necessary to map the plane z so that the unit circle is unchanged and the image of the 
secondary profile is mapped onto a quasi circle. This is done by 

,o 31 j zAl Li = ( z - z t> * - i/£f» y 

\ z - z Ns Z - l/z* Ns ) 

where r = n/(2n — tf s ), S s being the trailing edge angle of the secondary profile, zt s is the location on the 
plane z, and z^ s is a point inside the secondary profile in the vicinity of the leading edge. The constants 
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Fig. 2.7. Derivative of the circulation with respect to x s - 
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Fig. 2.8. Circulation as a function of x s and 9 S : the circulation has a local minimum. 


and are found as follows. Let 


(3.4) 


and 

(3.5) 


f(z) = 


f z Z T s z 1 / Zj'g \ 

\z-z Ns z- 1 /z* N J 


9(0 = 


£ - £ ~ 1 Ht 

i/c N 


The mappings f(z) and g(£) have the singularities df/dz = dg/d£ = 0 for f(zpj s ) = g((,iv) and f(zT S ) — 
#(£t)* These equations were solved iteratively in [6], whereas we found that they can be solved in a closed 
form. 

At last, by the Theodorsen-Garrik transform it is possible to map £ on two concentric circles on the 
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Fig. 2.9. Derivative of the circulation with respect to x 3 . 



plane z c : 

oo 

(3-6) ^ = z c exp{^[(-a' + ib'^Rz^ + (a' + ib'^R/z^}} 

j = 0 

where z c = J?cxp(i0). The radius R and the coefficients a' and b j are obtained by a trial and error process 
based on an FFT, as for Eq. (3.2). For example in Fig.( 3. 2-3. 3) the pressure isolines and the streamlines, 
relative to the flow around two wing sections computed by this method, is presented. 

3.1. Divided differences derivative. The set of transformations we described in the previous section 
is not easily differentiated explicitly, since the Theodorsen tranform and the Theodorsen- Garrick transform 
are defined implicitly, and the coefficients are determined by an iterative process. For this reason we decided 
to differentiate the mapping from the physical plane to the plane of the concentric circles by divided differ- 
ences. Then, we used the results of the analytic differentiation of Section 2 in order to compute the gradient 
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Fig. 2.11. Derivative of the circulation with respect to x s . 


dG dthsQ E 



1 0 ths 


Fig. 2.12. Derivative of the circulation with respect to 0 9 


of the circulation with respect to the three parameters which define the position of the flap with respect to 
the main airfoil in the physical space, i.e., (X S ,Y S ) the location of the leading edge of the flap and 0 its 
angle of attack. 

As we mentioned in Section 2, all of the variables (2.11) defining the solution on the plane of the adjacent 
circles with focal points in ±1 are functions of X s , Y s and 0. We computed the derivative of such functions 
by divided differences: for example we gave a small increment to X 8 and evaluate the increments of the 
variables (2.11) by means of the chain of mappings described in the section above. The differential is then 
obtained as the ratio between the increments of variables (2.11) and the increment in X s . 

As an example, we present the results obtained for the case of two adjacent NACA-0012 airfoils. The 
incidence of the flow at infinity is 10 degrees and the chord of the flap is 1/10 the chord of the main airfoil; 
also, the range for X s , Y s and 1 < X s < 1.3, -0.1 < Y s < 0.1 and 0 is and -5 < 0 < -2. These intervals 


Fig. 3.1. Plane of the two airfoils . 



are subdivided on a mesh of 10 points. As it is seen in Fig. 3.4, there exist a non-monotone region of the 
circulation manifold. The circulation isosurfaces are plotted for increasing values of lift from top to bottom 
(0.71 to 0.74). In Fig 3.5 the gradient is shown with respect to 0 computed as explained in the present 
section. The gradient values of the plotted isosurfaces range from 0 to 8 - 10“ 3 from left to right. 

The features of the circulation manifold are influenced by * he dimension of the flap as well as by the 
incidence of the flow at infinity. An abundance of different behaviors may be obtained changing these 
parameters, as in the two adjacent circles case. 

4. Conclusions. We presented two optimization test cases. The first concerns the computation of 
the circulation functional for two adjacent circles, when the relative position and the trailing edge of the 
secondary circle are varied. The gradient with respect to those variables was computed analytically and 
results have been presented for several geometrical configurations. 
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Fig. 3.3. Streamlines. 

10 



Fig. 3.4. Circulation isosurfaces: increasing values from top to bottom. D-x is the increment in X a , D-y in Y a and D-a 
in ©. 


The second test case was related to the computation of the circulation functional for two adjacent airfoils, 
as the secondary airfoil is displaced. Since the mappings involved to compute the theoretical solution are 
defined implicitly, we decided to compute the gradient of such transformations by divided differences. Results 
showing the nature of the circulation manifold were presented. 

In summary, it was shown that the space of solutions for such problems include non trivial situations of 
local minima and maxima. Such richness of behaviors and the relatively complicated geometrical optimiza- 
tion, which in the case of numerical solution involves remeshing and accurate definition of the gradient for 
moving boundaries, should be an interesting feature to test optimization algorithm of different nature. 

The codes developed for this study are available at ICASE - NASA Langley Research Center, Hampton, 
VA. 
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Fig. 3.5. Gradient isosurfaces. Gradient values from 0 to 8 • 10 3 from left to right . 

Aknowledgments. We would like to thank Manuel D. Salas, ICASE Director, for pointing out that 
the assumptions in [3] were erroneous, and hence, for motivating this work. Also, we are pleased to thank 
Professor Gino Moretti, N.Y. Polytechnic, for several valuable discussions resulting in original hints and 
remarks about the mapping of two profiles onto two circles and about an alternative derivation of the elliptic 
functions with respect to their period. 


REFERENCES 

[1] M. Lagally, Die reibungslose Stromung im Aussengebiet zweier Kreise , ZAMM, Band 9, Heft 4, Aug. 

1929, pp. 299-305. 

[2] A. Iollo AND M.D. Salas, Optimum transonic airfoils based on the Euler equations , ICASE Report 

No. 96-76, 1996. 

[3] B.R. Williams, An exact test case for the plane potential flow about two adjacent lifting airfoils , R.A.E. 

Report No. 3717, 1971. 

[4] C. Ferrari, Sulla trasformazione conforme di due cerchi in due profili alari , Memorie della Reale 

Accademia dclle Scienze di Torino, Serie II, Vol LXVII, 1930, pp. 1-15. 

[5] F. Tricomi, Funzioni ellittiche , Monografie di Matematica Applicata CNR, Zanichelli, 1951. 

[6] D.C. Ives, A modem look at conformal mapping including multiply connected regions , AIAA Journal, 

14, No. 8, 1976. 

[7] N. Marco, C. Godart, J.A. Desideri, B. Mantel, and J. Periaux, A genetic algorithm compared 

with a gradient-based method for the solution of an active control model problem , INRIA Research 
report RR-2948, INRIA, France. 


16 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services. Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Papenvork Reduction Project (0704-0188), Washington, DC 20503. 


1. AGENCY USE ONLYfCeave blank) 2. REPORT DATE 

July 1998 


REPORT TYPE AND DATES COVERED 

Contractor Report 


4. TITLE AND SUBTITLE 

On the circulation manifold for two adjacent lifting sections 

5. FUNDING NUMBERS 

C NAS1-19480 
C NAS1-97046 

6. AUTHOR(S) 

Luca Zannetti 
Angelo Iollo 

WU 505-90-52-01 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Institute for Computer Applications in Science and Engineering 
Mail Stop 403, NASA Langley Research Center 
Hampton, VA 23681-2199 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

ICASE Report No. 98-27 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, VA 23681-2199 

10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

NASA/CR- 1998-208447 
ICASE Report No. 98-27 

11. SUPPLEMENTARY NOTES 

Langley Technical Monitor: Dennis M. Bushnell 
Final Report 
Submitted to ZAMM 

12a. DISTRIBUTION/AVAILABILITY STATEMENT 

12b. DISTRIBUTION CODE 

Unclassified Unlimited 
Subject Category 34, 64 
Distribution: Nonstandard 
Availability: NASA-CASI (301)621-0390 



13. ABSTRACT (Maximum 200 words) 

The circulation functional relative to two adjacent lifting sections is studied for two cases. In the first case we 
consider two adjacent circles. The circulation is computed as a function of the displacement of the secondary circle 
along the axis joining the two centers and of the angle of attack of the secondary circle. The gradient of such 
functional is computed by deriving a set of elliptic functions with respect both to their argument and to their period. 
In the second case studied, we considered a wing- flap configuration. The circulation is computed by some implicit 
mappings, whose differentials with respect to the variation of the geometrical configuration in the physical space are 
found by divided differences. Configurations giving rise to local maxiira and minima in the circulation manifold are 
presented. 


14. SUBJECT TERMS 

potential flow; elliptic functions; optimal shape design 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


18. SECURITY CLASSIFICATIO 
OF THIS PAGE 

Unclassified 


15. NUMBER OF PAGES 

21 


16. PRICE CODE 

A03 


19. SECURITY CLASSIFICATION! 20. LIMITATION 
OF ABSTRACT I OF ABSTRACT 


551 


Standard Form 298 (Rev. 2 

Prescribed by ANSI Std. Z39-18 
298-102 




















